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Abstract: In this paper we analyse three models of the early universe, for which the 
respective mechanisms for generating the curvature perturbation are considered disparate. 
We find that in fact the mechanisms are very similar, and hence explain why they give 
rise to a large non-gaussianity. We show that the mechanism for generating the primordial 
curvature perturbation, and hence the observable non-gaussianity, is similar in both the 
Curvaton and Modulated Reheating models. In both cases the model can be written in 
terms of an energy transfer between the constituting fluids. We then show that this is also 
true for the mechanism of generating the curvature perturbation by symmetry breaking 
the end of inflation. We then relate this to the non-gaussian contribution to the curvature 
perturbation and find that it is inversely proportional to the efficiency with which the 
curvature perturbation is transferred between the fluids. For the first time, we generalise 
models of modulated reheating to allow for a non-linear energy transfer rate. 
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1. Introduction 

One of the greatest successes of inflationary cosmology is to provide a mechanism to gen- 
erate the spectrum for the primordial density perturbations, in excellent agreement with 
recent observations of the Cosmic Microwave Background (CMB) Q and galaxy surveys 
(e.g. 0]). In the most popular single field inflation models, the same field that drives in- 
flation is also responsible for the generation of the spectrum during inflation. However, in 
recent years models in which these tasks are separated have become prominent, in particu- 
lar the curvaton scenario || |], modulated reheating || [7], || H, and an inhomogeneous 
end of inflation [l^, |l3|, Q (here the separation is incomplete). Although these 

models were usually considered to be very different, taking the mechanism that generates 
the primordial spectrum of perturbations to categorise the models we show that these mod- 
els are in fact very similar, as the generation mechanism of the perturbations is the same. 
This also explains why all three models can give rise to large non-gaussianity. 

The same observations cited above also show that the spectrum of primordial density 
perturbations is predominantly gaussian in nature with a possible minor deviation JToj , |16[ . 
Single field canonical models of inflation are known to give a non-gaussian contribution 



which is slow roll suppressed [17]; this level of non-gaussianity is undetectable by current 
technology. It has been shown that if the field that generated the primordial curvature 
perturbation interacts with another fluid, such as in the curvaton (13, Il9| or modulated 
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reheating || [jj ^, or if inflation is driven by more than one field, see e.g. 
then the non-gaussianity can be enhanced to a detectable level. 

In this paper we evaluate the curvature perturbation £ and the non-linearity parameter 
/nl for models of modulated reheating, the curvaton model, and the inhomogeneous end 
of inflation. We investigate the relationship between these models and highlight the un- 
derlying mechanism of generation and the conditions for enhancing /nl- Although at first 
glance these models appear to be very disparate, we show that the generation mechanism 
in all three cases is similar: all three models can be written in terms of an energy transfer 
between the constituting fluids; the energy transfer itself is controlled by a light scalar field; 
the fluctuations of the scalar field are then "inherited" after its decay by the curvature per- 
turbation. In the curvaton and modulated reheating models, the curvature perturbation 
can be parametrised by the efficiency with which the curvature perturbation is translated 
from the field that generates the primordial curvature perturbation to the radiation field 
and hence is the determining factor in /nl- The same is true for the inhomogeneous end 
of inflation case, however the efficiency parameter here refers to the efficiency by which 
the curvature perturbation is translated from one scalar field to the one that generates the 
primordial density perturbation. 

For the first time we consider an energy transfer which depends non-linearly on the 
energy density in the modulated reheating scenario. We show how significantly more gen- 
eral scenarios lead to a curvature perturbation with the same functional form but modified 
coefficients, which affects the value of /nl but it still becomes large in the same limit as 
in the standard case of a linear transfer. Unlike previous calculations, we do not resort 
to solving the equations of modulated reheating numerically, which makes the derivation 
more transparent. 

This paper is organised as follows: in Section (§) we present the equations governing 
the energy transfer between fluids, in Section (||) we evaluate the curvature perturbation 
for models in which a scalar field decays into radiation, in Section (||) we calculate the 
curvature perturbation using a different approach to that derived in previous work for the 
inhomogeneous end of inflation. We present the results for the non-gaussianity parameter 
/nl in Section (||) and discuss them in Section (|6|). We use the convention of decomposing 
the curvature perturbation in terms of the first order £j and second order £2 contributions 
as follows 

C = Ci + |c 2 , (i.i) 

where Ci is the linear (gaussian) perturbation and C2 is a gaussian squared, see e.g. Ref. p3|| . 
2. Energy Transfer 

The line element for a Friedmann-Robertson- Walker spacetime including scalar perturba- 
tions, without yet specifying a particular gauge, is 

ds 2 = -(1 + 2<p)dt 2 + 2aB ti dtdx i + a 2 [(1 - 2ip)5 tj + 2E >ij ] dx l dx j , (2.1) 



20, 21, 22], 
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where we use the notation and conventions of Refs. [23| with the curvature perturbation, ifi, 
the lapse function, <ft, and scalar shear, a 2 E — aB, and where 8ij denotes the flat background 
metric and X, = dXjdx 1 . 

Allowing for the exchange of energy between fluids, the equations governing the evo- 



lution of the density perturbation are given by Refs. [24, p3|] . From the local energy- 
momentum equation, 

= Q U ' (2 - 2) 

we have the evolution equation for the energy density of a particular fluid in the background 

p a + 3H{ Pa + P a ) = Q a , (2.3) 

and at first order in the perturbations on superhorizon scales 

Sp a + 3H(5p a + SP a ) - 3(p a + P a )^ + Q a ^ -8Q a = 0, (2.4) 

where p is energy density, P is pressure, subscript a denotes the fluid species, Q a is the 
energy transfer to the a fluid and T^ v is the stress-energy tensor. 

Note that in Sections ( |Q| ) to fl3.3| ) we use the longitudinal gauge, a 2 E — aB = 0, and 
we ignore anisotropic stress so <p = tp, whereas in Sections ( |3.4| ) and (||) we use the flat 
gauge, ip = 0, for convenience. 

3. Models in which a scalar field decays into radiation 

We now consider three different forms for the energy transfer: in model ([O]) Q m is a 
linear function of p m , in model ( |3.2|) it is a higher order function of p m and in model 
Q3.3|) it is a function of a combination of powers. These three models are all examples of 
modulated reheating, in which the decay rate from the inflaton to radiation is modulated by 
a subdominant scalar. Since the Univere is matter dominated before decay and radiation 
dominated after decay, changes in the decay time lead to changes in the expansion and 
hence effect the curvature perturbation ||. In section ( |3.4|) we present the generic results 
for the curvaton scenario, in which the rate of decay is homogeneous. 

3.1 Q m as a linear function of p m 

We model the energy transfer Q a according to the standard choice 

Qa = -Tp a , (3.1) 

where V is the decay rate of the a field into radiation. In this section we consider modulated 
reheating in which the decay rate of a light scalar field ip varies, that is V = T(ip), and 



in Section (3.4) we analyse the curvaton (a) model in which T is taken to be constant. 
For both these models the background evolution equations for an oscillating scalar field 
(i.e. matter type) fluid and radiation are [^. Q 

-3Hp m - Tp 

Pr = -4:Hp r + Tp m , (3.2) 
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and we assume that the decay is into Standard Model radiation. Evaluating Eq. fl2.3|) for 
the oscillating scalar field, noting that P m = 5P m = and 5Q m = —5Tp m — T5p m , which 
follows from (|3.2| ), we get 

^ + 3H ^l-34> + r<i > + r^ + 5T = 0. (3.3) 

Pm Pm Pm 

Introducing the density contrast 5 a = Sp a /p a to rewrite the above equation in a more 
compact form, and using Eq. (EO) we get 8p m /p m = 5 m — 5 m (3H + T) and 



6 m = 3<f> - T(</> + 5r) , (3.4) 

where 5r = 6T/T. Similarly, for the radiation fluid, where P v = p r /3 and Q r = Tp m , we 
find 

6 r = 4<j> + T^(<f> + 5 T + 5 rn -5 I ). (3.5) 

Pr 

Finally we have from the time component of the Einstein equations G^ v = 8ttGT^ u on 
large scales [25] 

H<p + H 2 <f) = -^-( Pm 5 m + Pr 5 r ). (3.6) 

In Ref. H Modulated Reheating is studied in the so called 'purely forced' case, where 
the metric perturbations are sourced only by the perturbations in the rate of decay. To 
derive the curvature perturbation for this case we assume that the perturbations in the 
matter and radiation fields are subdominant, hence leaving us with two equations 

H(p mo d + H 2 (f> mod = , 

30mod = r(</> mo d + 8r) , (3.7) 
where the subscript "mod" refers to 'modulated reheating'. The above equations then give 

^mod = "f^FJ^r = -/3— . (3.8) 

Here f3 is proportional to the decay rate with /3 — 1 corresponding to instant reheating. 
The curvature perturbation on uniform density hypersurfaces is defined as 

Cmod = — f^mod ~ H— . (3.9) 
P 



where Sp/p = Yli^Pi/ Y2j Pj- Evaluating the RHS of Eq. ( |3.9| ) in the "forced case" 
(i.e. when 5p/p -C (ft) gives 

Cmod = — 0mod j (3.10) 

and finally, substituting Eq. (pH) into Eq. ( |3.10 ) we get 

Cmod = $Y ■ ( 3 - n ) 
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Assuming now that T depends only on one light scalar field, ip, (Ref. || considers a 
more general case) a simple Taylor expansion results in 



dT 1 d 2 T 

r = r + — 5^ + -—^, (3.12) 

dip* 2 dpi 

where the subscript * refers to horizon exit, this gives the curvature perturbation as 

Cmod = ^ (r v j^ + , (3.13) 

where and T Vfipf are the first and second derivatives with respect to </?*. 

In sections ( |3.4| ) and (Q) we define the efficiency as the ratio of the first order curvature 
perturbation of a field prior to its decay with respect to its value post decay, and hence we 
construct similar parameters for this model for completeness. The field p* can be related 



to a curvature perturbation of = (<£>'* + Vtp t )5p*/(3ip\ 2 ) = b5p*, see e.g. Ref. [26], then 
the efficiency parameter is 

Cmod 1 @^,<p* /, -1 i\ 

Cmod 2 /5r,iz>*c5» /„ -. r x 

C2 = -g- = ■ (3-15) 

The overall curvature perturbation for this model can now be written as 

Cmod = bcxSip* + b 2 c 2 5(fl . (3.16) 
3.2 Q m as a higher order function of p m 

In this case we define the energy transfer as a function of a higher power of the energy 
density, similar to the model recently considered in the context of dark energy decay in 

Qm = Tp n m . (3.17) 

In order to avoid singular behaviour in the limit of small p m we require n > 0. 

Following the same steps as in section ( |3.1[ ), i.e. only considering bp = 0, we find that 
the curvature perturbation on superhorizon scales has the same form as Eq. (3.13), with /3 
now also a function of the energy density: 

To 11 - 1 

ZH + Tp^ 1 1 ; 

In order for p m to decay then, the decay rate must exceed the expansion rate of the Universe, 
Pm *C 10 -12 in Planck Units, and since H decays as pll" 2 , then in order for H/Tp 1 ^ 1 < 1 
we require n < 3/2. 
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3.3 Q m as a function of a combination of powers of p m 

In this case we define the energy transfer function as 

Q m = r (ApP m + BfflJ , (3.19) 
where < p < q. Then C has the same functional form as Eq. ( 3.13j ) but with (3 redefined 



as 



r (a fa 1 + BpV 

P = '—^r • (3-20) 

3H + r [A fa 1 + BpV 



3.4 The curvaton scenario: 5T = 

We now consider the second scenario, where 5T = 0. This is the case in the curvaton 
model. The energy density of the curvaton is highly subdominant compared to the inflaton 
field during inflation, but afterwards, while it oscillates about the (quadratic) minimum of 
its potential, its energy density decays like matter. This energy component grows relative 
to the radiation density generated by the inflaton decay products until the curvaton also 
decays into radiation [§, f|, [j. In this scenario the curvature perturbation is solely generated 
from perturbations in the curvaton field. Using a quadratic potential for the curvaton, the 
curvature perturbation in this model was found to be 

\ 2 \ 

(3.21) 

where a has been split as a = <7o* + #01*, and r\ parametrises the contribution of the a 
field to the overall curvature perturbation post a decay 

Ccurv l,post a decay /„ rir\\ 

7*1 — — . yo.Az) 

S,curv l,pre a decay 




Using second order perturbation theory Ccurv is given by []19| 

Ccurv = ~ri + - ri + -r 2 , (3.23) 

3 do* o \ 3 y V <7o* 



where 

Ccurv 2, post a decay 

c 5- 

^cur\ 

4. Inhomogeneous End of Inflation 



^curv z,posi cr aecay / ^ a\ 

r 2 = — 2 * y - . (3.24) 

'curv l,pre a decay 



In previous work Jl0| , 11, it has been shown that the curvature perturbation can be 
generated at the end of inflation if there is an ultra light scalar field sub-dominant to the 
inflato during inflation. This second field does not play a role in the inflationary dynamics, 
but serves to perturb the inflaton trajectory from a straight line, potentially resulting in 
the generation of non-gaussianity. This concept was generalised to a two-field hybrid model 
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of inflation in [12| and was found to give a measurable /nl- The potential describing this 
model is given by: 



(f^ +gip a + ha 2 )^- + ^ X 2 , 



(4.1) 



where % ls the waterfall field which is held at zero during inflation, and /, g and h define 
the coupling of the inflaton fields to the waterfall field. In this case / 7^ g 7^ h and hence 
the surface at the end of inflation is best defined by an ellipse. For simplicity we assume 
an instantaneous decay of the scalar fields into radiation at the end of inflation, for a more 
detailed discussion of this point see Refs. [ 



In contrast to previous work carried out on this and similar models [12, [L0|, 11, 13, 14] 
which used the 5N formalism, we use the perturbative approach in calculating /nl- Starting 
with the definition of the curvature perturbation at horizon exit, 



Co 



-H 



5p c 



where the subscript "e" refers to the end of inflation. We take the a 
split ip e as (p e = ipo e + 5ip e , and thus, 



(4.2) 

trajectory and 



2% \ <A)e V<A)e. 



(4.3) 



In order to evaluate the spectrum and bi-spectrum for this model we will need to 
calculate terms of the form < 5^p c (ki)5ip c (k.2) >, but since we only know their values at 
horizon exit we re-write 5(p e in terms of S(p*: 



op* 2 



(4.4) 



where <// e = 8p c /da = -g/2f, <p c = -(2h - g 2 /(2f))/(2f<p c ), dip e /dip* = e^ N @, rj s is 
the second derivative slow roll parameter and N is the logarithmic ratio of the scale factor 
at the end of inflation with respect to its value when scales of cosmological interest exited 



the horizon [29]. 

Assuming that the curvature perturbation is generated predominantly from the a field, 
then Eq. (|4.3j) becomes 



5a* + 



1 



'2 " \ r 2 



(4.5) 



and since we are mainly interested in highlighting the underlying physical mechanism, we 



have used a simplifying [18] calculation to get the second order contribution, which explains 



why our results do not match those of Ref. [12] at second order exactly. 
We now define the parameters 



(a 



Cip 1 



(4.6) 
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and 



i Co- 2 Tj s _2n,JV ( 12 , ' 

V 1 Z V 



(4.7) 



in analogy with the parameters c\ and C2 above. Therefore the curvature perturbations 
generated by the a field can be written as 




4-—( — * 



(4.8) 



This model can be recast in the form of two interacting fluids by writing Eqs. (| 
and ( |2.4D in the slow roll approximation for this case 



-Qa = 
-5Q a 



-V a & = . 
V 2 



(4.9) 



showing more clearly the relation to the previous models and hence explaining the similar 
results in terms of the non-linearity parameter, as detailed in the next section. Note that 
we are only modelling the energy transfer between fields and not the decay of the field into 
radiation. We assume an instantaneous decay into radiation at the end of inflation, for a 

13 1 . The energy of the (p field is transferred to the 



discussion of this point see Refs. 



radiation at the end of inflation and the perturbation of this energy transfer is dominantly 
controlled by the perturbation of the a field. After the decay of the ip field, the energy of 
the radiation was endowed from the (p field but the perturbation is inherited from the a 
field. 



5. The Non-Gaussianity Parameter /] 



NL 



From Eq. ( ^16|) we find that the non-linearity parameter /nl for the modulated reheating 
case is given by 

, 5 c 2 
= 6 ' 



which in terms of j3 is 



p 5 r r (yj,)/)* 

/NL = T^~B~ 



12/3 



Using the functional form V = T\{\ + Sip^/ip*) 2 , following ||, Eq. (^2) gives 



(5.1) 
(5.2) 



/: 



NL 



24/3 



(5.3) 



and /3 depends on the decay rate. This inverse dependence on j3 is also valid for the new 



models discussed in Sections (3.2) and (|3.3| ). 

The non-linearity parameter /nl for the curvaton model is derived in the simplest case 



from Eq. JOlp to be M 



NL 



4ri 



(5.4) 
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and using second order perturbation theory |T{| 



5 5r 2 5 



(5.5) 



4ri 6rf 3 



Eqs. ( |5.4| ) and ( [5.5D were evaluated for the potential V(<r) = m 2 a 2 /2. 
Finally we have for the inhomogeneous end of inflation: 

5 6 2 



(5.6) 




6. Discussion 

In this paper we have studied models of the early universe dominated by a scalar field that 
then decays into radiation, and a model in which two scalar fields drive inflation where 
the curvature perturbation is generated primarily at the end of inflation. Using the energy 
transfer Q = Tp to characterise the models, the former case can be further split into models 
with a homogeneous T and varying p, Curvaton models, and models with a varying decay 
rate F and homogeneous energy density p, Modulated reheating. 

We then evaluated /nl for the simplest case of a homogeneous T, and for various 
functional forms of the energy transfer parameter Q v . We found that in the varying F case 
/nl is inversely proportional to the decay efficiency /3, even in the new case we consider 
of a non-linear energy transfer, Q = Tp n or Q = F (Ap?n + Bptn), which parametrises the 
rate at which the scalar field decays into radiation, as well as the energy density of the 
scalar field prior to its' decay. Reparametrising our equations shows that /nl is inversely 
proportional to cf , where c\ is the ratio of ( prior to the decay of the light scalar ip to it's 
value post (p decay. Similarly in the case where the curvature perturbation is generated 
by the inhomogeneous end of inflation, we found that /nl is inversely proportional to the 
square of the parameter b\, the ratio of the first order curvature perturbation sourced by a 
to the curvature perturbation sourced by ip at horizon exit. Previous work [TIJ] (which we 
reproduced here) already showed a similar result for the curvaton model; /nl is inversely 
proportional to r%, the ratio of the first order curvature perturbation after the decay of a 
with respect to the value of the curvature perturbation prior to its decay, Eq. ( ^22|) . 

Our results show that the mechanism of generating the curvature perturbation in these 
models is similar, and can be modelled as energy transfer between interacting fluids: the 
fluctuations are generated in the first fluid, a massless scalar field, and are then inherited 
by the second one into which it decays. This therefore explains why the condition for 
enhancing /nl is similar; it is the efficiency with which the energy density or curvature 
perturbation is translated between the fluids that determines the level of non-gaussianity. 

Although the three models we have discussed are physically very different and neither 
share the same fields nor generate the curvature perturbation during the same epoch, in 
all cases the curvature perturbation is related to an energy transfer. This energy transfer 
is from a light scalar field into some other fluid into which it (and any other fields present) 
decay, see Table (Q). We have explicitly shown that in all three cases it is the efficiency 
with which the energy density or curvature perturbation is translated between the fluids 
that determines the level of non-gaussianity. 
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Table 1: We summarise the generation of large non-Gaussianity in the three models analysed in 
this paper. The second and third column show the transfer of energy and the fourth column shows 
the source fluid which modulates the energy transfer and generates the curvature perturbation. The 
fifth column shows the source of the dominant background energy density. Large non-Gaussianity is 
only possible when the source of the background energy and that of the perturbations are different. 
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